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Recently the Lattice Strong Dynamics Collaboration reported a reduction of the electroweak S 

parameter for SU(3) gauge theory with Nf = 6 fermions in the fundamental representation, com- 

pared to scaled-up QCD. Here I provide additional details of our calculation. I discuss our use 

of conserved lattice currents; the relation to vector-axial parity doubling; finite- volume effects; 

and the sensitivity of our results to the number of fermion doublets with chiral electroweak cou- 

plings. Results presented here include additional data, and do not affect our previously -published 

^ conclusions. 
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1. Introduction 

The application of lattice gauge theory to strongly-interacting physics beyond QCD is at 
present a very active field [1]. While much of the current interest is motivated by the possibil- 
ity that new strong dynamics may play a role in electroweak symmetry breaking [2, 3], improving 
our general understanding of strong dynamics is an important theoretical goal in its own right. 

The standard picture of strongly-interacting Sl](N) gauge theories is that as we increase the 
number Nf of fermions in a given representation, an infrared fixed point will develop at some 

(c) (c) 

critical Nj- ^ . For Nf > Ny (up to the loss of asymptotic freedom) the system is IR-conformal. 

Approximately-conformal systems with Nf < A/j^^^ may possess the dynamical scale separation that 
characterizes "walking" theories, as well as parity doubling between vector (V) and axial-vector 
(A) spectra that can reduce the electroweak S parameter to phenomenologically viable values [4]. 

The Lattice Strong Dynamics Collaboration approaches these questions by using QCD as 
a baseline. We consider SU(3) gauge theory and steadily increase the number of fundamental 
fermions, comparing our results against the familiar case Nf = 2. We use computationally ex- 
pensive domain wall fermions for better control over lattice artifacts. Our first studies focused on 
the Nf = 6 model, which while not truly walking exhibits some of the associated phenomena: by 
matching IR scales between Nf = 2 and Nf = 6 calculations, we observed an enhancement in the 
Nf ^ 6 chiral condensate [5] and a reduction of the S parameter relative to scaled-up QCD [6]. 
Here I provide additional details of our S parameter calculation that were not discussed in Ref. [6]. 
Results presented here also include additional data, and do not affect the conclusions of Ref. [6]. 

We can identify three main ingredients in our expression for the S parameter. 



The term A5w accounts for the three Nambu-Goldstone bosons (NGBs) eaten by the and Z, 
and is discussed in detail by Ref. [6]. In Section 2 I review our calculation of the transverse V-A 
polarization function Ily-A{Q^)^ and relate it to the vector and axial spectra in Section 3. Finally, 
Nd is the number of doublets with chiral electroweak couplings; in Section 4 I show how it affects 
our results for the S parameter. 

2. Currents and correlators 

On the lattice, the transverse V-A polarization function Ily-A{Q^) is determined from 



Here Q = 2sin(7rn/L) are lattice momenta, while Q = 2Kn/L\ these are spacelike = —q^ > 0. 
The current correlators mix two types of domain wall currents. V^^ and A^^ are non-conserved 



S^AtiNd lim 



ny_A(Q')-A5w. 



(1.1) 




(2.1) 
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"local" currents defined on the domain walls; in terms of five-dimensional fermion fields ^{x,s). 



1 rr^ 



V^'^ix) = - {^(;c,Lj - 1)/'(1 + ^)^^>{x,Ls - 1) +»I'(x,0)/'(l - ^)^'V{x,0)] 
A^"{x) = ^ {^(x,L, - \)f{\ + ^)^m{x,L, - 1) -»I'(x,0)/'(l - ^)Tf^'{x,0)] 



(2.2) 



The conserved currents Y^" and s^^'^ are point-split, and summed over the fifth dimension: 



r^-{x)='^Y.f-{x,s) £\ign/',_:^ (2.3) 



(x, 5) = - {^(x + M , 5) ( 1 + 7^ )f/,V 5) - ^(x, 5) ( 1 - / ) f/,,^ T«^(x + M , ^) } . (2.4) 

The Fourier transform in Eqn. 2. 1 involves (x + /i /2) because the conserved currents are point-split 
on the link + /i). The flavor matrices are normalized to Tr [t^t^] = 5^^/2. 

Although the conserved and local currents must agree in the continuum limit, at finite lattice 
spacing only the former satisfy a Ward identity {Q^Ylyy = 0, Fig. 1). Because the correlators 
involve both currents, Eqn. 2.1 includes the renormalization factor Z, which we compute non- 
perturbatively, Z = 0.85 (0.73) for Nf = 2 (6). Our chiral lattice fermions ensure that Z^Za^ Zy. 
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Figure 1: On every configuration, Q/^Ilyy = when one conserved current is used in each correlator (left), 
but not when only non-conserved local currents are used (right). The horizontal offsets around each value 
distinguish different v. 



In principle, it would be best to work entirely with the conserved currents and 
instead of using the mixed correlators in Eqn. 2.1. In practice, evaluating conserved-conserved 
correlators such as (y^^(x)y^^(O)) requires ^{Ls) inversions, increasing the computational cost 
of the calculation by roughly an order of magnitude. As emphasized in Ref. [7], lattice artifacts 
cancel in the V-A difference of the mixed correlators, allowing us to use these less expensive 
quantities. This is illustrated in the left panel of Fig. 2: even though YL^^Qy ^ since V^^ and A^*^ 
are not conserved, [n^y (2^) — n^(2^)] 2v ^ 0. In the right panel, we see that this does not hold 
if we use only local currents in the correlators. 
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Figure 2: On every configuration, lattice artifacts Yl^^Qy ^ cancel in the V-A difference when one con- 
served current is used in each correlator (left), but not when only non-conserved local currents are used 
(right). The horizontal offsets around each value distinguish different \i. 



3. Parity doubling and finite volume effects 

Because chiral perturbation theory cannot reliably be applied to our A^/ = 6 calculations [8], 
we extract the slope ny_^(0) by fitting our data to a simple four-parameter rational function, 

This "Pade(l,2)" functional form has the correct asymptotic behavior ny_A(2^) ^ at large 
Q^, and also resembles the single-pole dominance approximation to the Y-A dispersion relation 



Rv{s)~Ra{s) 

5 + 



(3.2) 



{Fp is the pseudoscalar decay constant.) That is, with the single-pole dominance approximation 
R{s) — 12k^F^5{s — M^), this dispersion relation becomes 

f-f{pole) y _ Tp2 I Q^^V Q^^A /o ON 

nv_. (Q )--/> + - j^,, (3.3) 

which reproduces the form of Eqn. 3.1 when we apply the corresponding approximation to the first 
Weinberg sum rule, Fp ^ Fy —F^. Because the lattice data contain information about the entire 
spectrum, the fit parameters in Eqn. 3.1 do not directly correspond to the combinations of meson 
masses and decay constants predicted by the pole-dominance Eqn. 3.3. 

Uncorrelated fits of our data to Eqn. 3.1 produce stable results with jdof <C 1 as we vary 
the fit range. Our results for ny_^(0) are shown as colored points in the left panel of Fig. 3. 
The black points in that plot are pole-dominance predictions based on Eqn. 3.3. Both the direct 
fit results and the pole-dominance predictions show a reduction for Nf ^6 compared to A// = 2 
at light pseudoscalar masses Mp < Myo, where Myo is the vector meson mass in the chiral limit. 
The pole-dominance predictions are systematically lower than the direct results, consistent with 
the expectation that states neglected by the single-pole dominance approximation would provide 
additional positive contributions. 



4 



S parameter and parity doubling below the conformal window 



David Schaich 




Figure 3: The slope of ny_A(G^) at = 0, plotted versus mI/M^q. Left: Nf = 2 and 6 results on 32^ x64 
volumes from direct fits to Eqn. 3.1 (colored), compared to pole-dominance predictions (black). Right: 
Nf = 6 results on 16^ x 32 and 32^ x 64 volumes. 

The lightest Nf = 2 points in Fig. 3 are empty because they correspond to a fermion mass m 
so small that finite-volume effects may be significant. Finite- volume effects are a concern for the 
S parameter calculation because they can produce spurious parity doubling that artificially reduces 
Hy_j^{0). This is illustrated in the right panel of Fig. 3 for A// = 6 calculations on 16^x32 volumes: 
ny_A(0) as m — > 0, which would naively suggest a negative S parameter from Eqn. 1.1. The 
associated distortion of the spectrum provides clear evidence that this is merely a finite-volume 
effect: as m decreases, the 16^ x 32 pseudoscalar mass Mp freezes around Mp ^ 1.2MyQ, which is 
not the case for the 32^ x 64 results also shown in the plot. 

Returning to the lightest Nf = 2 points, the pole-dominance prediction for lly_^{0) decreases 
due to spurious parity doubling from finite- volume effects. However, we do not see a similar re- 
duction in the direct fit result. Instead, this point clearly continues the trend established at heavier 
masses, and the corresponding Nf = 2 results for S (Fig. 4, below) reproduce the prediction ob- 
tained by scaling up QCD phenomenology, lim^2^Q5 = 0.32(3) [4]. This suggests that the Fade 
fits may be less sensitive than spectral quantities to these finite-volume effects, increasing our con- 
fidence that the reduction observed for A// = 6 is physical. 

4. 5 parameter results 

Realistic models of dynamical electroweak symmetry breaking must produce exactly three 
massless NGBs to be eaten by the and Z. Any additional pseudo-Nambu-Goldstone bosons 
(PNGBs) must acquire masses from standard-model and other (e.g., extended-technicolor) interac- 
tions in order to satisfy experimental constraints. On the lattice, however, we perform calculations 
with Nf — l degenerate massive PNGBs. When we use Eqn. 1.1 to determine the S parameter from 
the lly_^{0) results shown in Fig. 3, the A^^^ term removes the contribution only of the three 
would-be NGBs. (To be more precise, the /s = NGB does not contribute, and ASsm cancels 
the contribution of the I/3I = 1 pair.) The remaining Nj — 4 PNGBs introduce chiral-log terms 
oc log [Myo/Mp] that would diverge in the chiral limit Mp 0. 

Fig. 4 presents our S parameter results for Nf = 2 and 6, considering two possible values of 
Nd for Nf = 6. The plot on the left presents the case in which every fermion possesses chiral 
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electroweak couplings, No = Nf/2 = 3. The minimal case in which only a single doublet has 
chiral couplings (No = 1) is shown on the right. In both cases the A// = 6 results show a reduction 
compared to rescaling Nf = 2, before diverging in the chiral limit. With No = I the S parameter 
can be significantly closer to the experimental value S ^ —0. 15(10) for M^^^^ ^ 1 TeV [9]. 
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A Nf=2 

mNf=6 



0.0 as 



Figure 4: S parameter for Nf = 2 and 6, for the maximum No = 3 (left) and minimum No = I (right). The 
bands correspond to fits explained in the text. 

To guide the eye, we include in Fig. 4 simple linear fits accounting for the A^/) -dependent 
chiral-log divergence that remains for Nf > 2. We fit the lightest three solid points to the form 

1 



5 = A + 5x + 



127r 



log(iA) 



where x = Mj /My^ and tJ counts the pairs of PNGBs with I3 7^ 0, 

fovNo^Nf/l f| = 27V/-3 forAt> = l. 



(4.1) 



(4.2) 



The blue Nf = 6 curves allow us to estimate the fermion mass m at which we could directly observe 
chiral log effects. The necessary m is too small for us to explore on our present 32^ x 64 volumes. 
Again, in a realistic phenomenological context, we must have only three massless NGBs, with 



Nf—4 massive PNGBs. To estimate a definite value for the Nf = 6S parameter in this situation, we 

^2 - 0.38M^Q 



can imagine freezing the masses of all Nj 



4 PNGBs at some finite value (such as Mp 
at the minimum of the No ^ I blue curve in Fig. 4), and then taking only the three NGBs to the 
0. A qualitative picture of this scenario is sketched in Fig. 5. 



chiral limit Mp 
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Figure 5: S parameter for Nf = 2 and 6 with No = I, imagining that we freeze the masses of all Nj —4 
PNGBs at Mj = 0.38M^q, as described in the text. 
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